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Method And Apparatus For Extracting Cerebral Ventricular System From 

Images 

Field of invention 

The present invention relates to a method and apparatus for extracting cerebral 
ventricular system information from images of one or more cerebral ventricular 
regions. In particular, the present invention relates to extracting the human cerebral 
ventricular system from radiological images. 

Background of the invention 

The human cerebral ventricular system contains four intercommunicating chambers, 
the left lateral ventricle, the right lateral ventricle, the third ventricle, and the fourth 
ventricle. The ventricles contain cerebrospinal fluid (CSF) and changes in CSF 
volume and shape are associated with several neurological diseases such as congenital 
anomalies, post-traumatic disorders, pseudo-tumors, neuro-degenerative diseases, 
inflammatory diseases, headache, cognitive dysfunction/psychiatric diseases and post- 
operative changes. Quantification of the degree of dilatation of ventricles is important 
in the diagnosis of various diseases and for measuring the response to treatment. It 
may also be important in predicting the prognosis of the disease process. Therefore, 
extraction and quantification of the ventricular system from medical images is of 
primary importance. 

Normally all of the ventricles are interconnected to enable circulation of cerebrospinal 
fluid. Hence, in an ideal case, a single region growing operation should be able to 
extract a complete ventricular system. In reality, however, this does not work and 
many factors should be taken into account such as noise, spatial and contrast 
resolutions of the scan. In particular, the method has to handle a partial volume effect 
causing "leakages" from the ventricles to the extraventricular space, and 
intraventricular and extraventricular pathology distorting the normal anatomy of the 
ventricular system. 

Conventionally, the standard methods employed for investigating the human cerebral 
ventricular system have been ventriculography and pneumoencephalography, each of 
which has relatively high morbidity and mortality With the advent of high speed 
Computed Tomography (CT) and three-dimensional (3D) magnetic resonance (MR) 
imaging (MRI), robust techniques are required to extract and quantify the ventricular 
system. 

Due to the great importance of the ventricular system, its extraction has been the 
subject of extensive research work resulting in methods such as region-growing 
assisted by morphological operations, thresholding, template matching, atlas warping, 
level sets, active models, and knowledge based methods. 

Schnack HG, Hulshoff PHE, Baare WFC, Viergever MA, Kahn RS, "Automatic 
segmentation of the ventricular system from MR images of the human brain," 
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Neurolmage 2001 vol. 14, pp. 95-104 describes an automatic algorithm to extract the 
brSf fr ,° m T1 " wei ^ 3-D-FFE MR images of the 

on^t™ fT r ^ UP ° n region 'S rowin g ^ mathematical morphology 
operators. One of the limitations of this algorithm is its dependence on th? coarse 
binary brain segmentation image and the coarse white matter segmentation image 
Although precise white matter segmentation is not necessary precise bra£ 

E^T ^i° n * iS t0 i ventricular s y stem extraction. Another limitation is that it 

is not able to avoid inclusion of non-ventricular CSF regions although some 

SST£ ^T^f f &e SyStem haS been ^corporated tto Se 

Fu^tLnr^ , i OCa ?° n ■? *" SCed point < s > 311(1 searc ^S subregions. 
Furthermore the algonthm descnbed is not able to extract the complete ventricular 
system as the fourth ventricle is not extracted. Several artificial bounSeT^ 

ediwtlH e ^ ext S ctmg &e *** ventricle resnl ^ * some TZZ 

° ° btam preci !f results - 111 ^ra°A> the algorithm is time comuming 
to run (about 5 to 20 mm on a Pentium3 450 MHz PC for a brain volume). 

Worth AJ, Makris N, Patti MR, Goodman JM, Hoge EA, Caviness VS, Kennedy DN 
Trecise segmentation of the lateral ventricles and caudate nucleus m^bnrin 
ZSLTSIf a TT* ay « ^ hist °8 rams '" IEEE Transactions on Medici 

fte E 1 V°- 2 ,' PP - 3 ° 3 " 310 d6SCribeS 311 automated me *od *> extract 

™™H tentacles and caudate nucleus on Tl-weighted coronal MR brain images of 
normal subjects. The algonthm derives, automatically, intensity thresholds from 
anatomical mformation and a histogram, and then extracts the outlines of tht laS 

Zwl'l C °? f° 1 nUCl6U f by generating contours. The algorithm uses 

toowledge guided rules and methods for locating automatically certain anatomical 
landmarks (for example, the corpus callosum), and employs the local threTholTfor 

Zl^ ?*? ^ md CaUdatG nucleus ' ™ e ^ orith ^ * able t exlS 

on y the lateral ventricles m normal subjects. It needs manual editing of the resulting 

vo W 2? ? OT ^ larities * automatically generated outlines caufed by the pS 
volume effect and low contrast. It takes about 50 minutes to extract the lateral 
ventncles and caudate for a brain volume data. 

Kaus MR Warfield SK, Nabavi A, Black PM Jolesz FA Kikinis R, "Automated 

SffSfSSi* ^ ^ ° f brain tUm ° rS/ ' Radi ° logy 2001 ' voL 218, noTpP 
586-591 discloses an adaptive template-moderated classification method to extract 

ES^ ^ tU ^ 0rS> ? 6 iterates statistical classification to a^gn 

ilent f d r n "i mear Ration to align a digital anatomic atlas to the 
patient data Statistical classification is used to divide an image into different tissue 
classes on the basis of signal intensity values. If different tisfues have the same or 
overlapping gray-value distributions, additional information about the spatial location 
of anatomic stmctures is derived from a registered anatomic atlas As sucn, me 
algonthm is dependent on the digital atlas and is time consuming to run (5-10 
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Bikinis R, Shenton ME Losifescu DV. McCarley RW, Saiviroonpom P, Hokama 
HH Robation A, Metcalf D, Wible CG, Portas CM, Donnino RM, Jolesz FA "A 
digital brain atlas for surgical planning model-driven segmentation and teaching" 
^ Jrf 90 ^ ° n VisuaUzati °n and Computer Graphics 1996, vol. 2, no. 3, pp 
232-241 describes an automated model-based segmentation system that uses a digital 
brain atlas to extract brain structure. When MR images do not provide sufficient 
contrast between various structures, a combination of automated and supervised 
segmentation methods along with the atlas information is used to define regions of 
interest (for example, the ventricular system) and to guide a segmentation algorithm 

f «°f^ e ^ 0anat0miCal knowled g e - The anatomical atlas is registered with 
classified 3D images and an elastic matching method is used for the projection of atlas 
information onto the data sets of subjects and patients. The resulting segmentation of 
the ventricular system depends on the accuracy of the digital brain atlas, and the 
accuracy of the matching method. The algorithm works best with images of normal 
subjects and is time consuming to run. 

Baillard C, Hellier P, Barillot C, "Segmentation of 3D brain structures using level sets 
and dense registration," IEEE Workshop on Mathematical Methods on Biomedical 
Image Analysis (MMBIA 2000), pp. 94-101 describes a co-operative strategy for the 
segmentation of parts of ventricles from 3D brain MRI images which integrates 3D 
segmentation and 3D registration methods. The segmentation is based on thelevel set 
formalism. Using an automatic registration method to initialise the ventricular 
stoicture surface, the algorithm iteratively propagates the closed 3D surface towards 
toe desired boundanes through the evolution of a 4D implicit function. The major 
hmitations of this method are that it can only extract parts of ventricles and that it is 
tune consuming. 

Cootes T, Taylor C, Cooper DH, Graham J, "Active shape models - their training and 
application, Computer Vision and Image Understanding 1995, vol. 61, no.l, pp 38- 
59 discloses a statistical technique for building deformable shape templates and uses 
these models to extract various organs, including the ventricles, from 2D and 3D 
medical images The statistical formulation provides global shape constraints and 
allows the model to deform only in ways implied by the training set. The shape of the 
models represent objects by sets of landmark points that are placed in the same way 
on an object boundary in each input image. The points can be connected to form a 
deformable contour. By exaniining the statistics of framing sets of manually-labelled 
medical images and using principal component analysis, a shape model is derived 
that describes the average positions and the major models of variation of the object 
points. New shapes are generated using the mean shape and a weighted sum of the 
major modes of variation. Object boundaries are then extracted using this "point 
distribution model ' by examining a region around each model point to calculate the 
displacement required to move it towards the boundary. These displacements are then 
used to update the shape parameter weights. There are several problems with this 
approach The technique is sensitive to the initial position of the atlas, that is, if the 
initial rigid alignment is off by too much, the elastic match may perform poorly The 
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presence of neighbouring features may also cause matching problems, for example 
the atlas may warp to an incorrect boundary. Without user interaction, the atlas can 
have difficulty converging to complicated object boundaries. 

Wang Y, Staib LH, "Boundary finding with correspondence using statistical shape 

Tot ' £ C o ee ^ 8 ffiEE conference of computer vision and pattern recognition 
iyy», pp. 338-345 proposes an approach for finding the boundary where the 
correspondence of a subset of boundary points using local shape features of a model is 
simultaneously determined. Statistical point models are constructed with shape and 
shape variation generated from sets of examples using principal component analysis 
of the covanance matrix. The model is then analysed in a Bayesian scheme to find 
shape parameters and pose parameters. The algorithm is used to locate the boundary 
of the lateral ventricles. Although the statistical point model is able to capture 
considerable variability for the lateral ventricular boundary, it is not easy to build the 
statistical model specific to all classes of the ventricular system. Furthermore this 
approach is tune consuming. 

In 1998, as described in http://www- s i g .enst.fr/tsi/ g roi ros ^n/active a way of 
u Utm ? ? 6 C , 0 *f P° ndence between 311 MR volume and the atlas was proposed in 
which structural information (as flexible spatial constraints) was taken into account 
The spatial constraints were formulated by the adoption of fuzzy set theory and 
information fusion theory. Segmentation approaches were not used globally but 
conditionally to regions of interest with imprecise limits. The calculation for 
correspondence between the MRI data and the atlas inferred a discrete deformation 
field constrained by object surfaces. Those recognised objects were relied upon to 
extract the whole ventricular system. This appears to be the first method (until 1998) 
to include the fourth ventricle. The aforementioned method needs to be validated on a 
large set of images and there does not appear to be any published information 
concerning its ability to distinguish the ventricles which would enable the method to 
be assessed. From the results presented, the algorithm does not appear to be able to 
aut^ t ati2i e iy VentriCUlar aCCUratdy ' ° r to maint ain the connectivity 

Sonka M, Tadikonda SK, Collins SM, "Knowledge-based interpretation of MR brain 
images, IEEE Transactions on Medical Imaging 1996, vol. 15, no. 4, pp 443-452 
describes a fully automated segmentation method to extract brain structures including 
the s ventricular system from MR images. The algorithm is based on a hypothesize- 
and-venfy principle and uses a genetic algorithm (GA) optimisation technique to 
generate and evaluate image interpretation hypotheses in a feedback loop The 
algorithm is teamed on 20 out of 28 MR brain images, with observer-defined contours 
of structures being used as prior knowledge and incorporated in the genetic algorithm 
(GA)-based image interpretation method. The method is tested on the remaining eight 
brain images, and can produce accurate labelling results of neuroanatomies 
structures. One limitation is that it can only handle brain with no gross anatomic 
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2SZJ?*J!f- add «° n ' manUal identification of contours of neuroanatomies 
structures is tedious, time-consuming and vulnerable to inter-personal variations. 

™ SClmable HiUDLG ' "Q^^g smaU changes in brain ventricular 
SSSJSS 8 T ngld re ? Stration " ^CAI 2001, pp. 49-56 describes the 
implementation of a non-ngid registration algorithm based on optimising normalised 
mutual mfonnation to extract the lateral ventricles from MR images, if e 
uses a free-form deformation (FFD) to model local deformation, with the FFD S 
constructed from a 3D tensor product of B-splines, and deformation being achieved 
by translating control points in steps along the direction with maximum gradient until 
either the magnitude of the gradient is less than or equal to a threshold, or a^e- 

n^ttrS ° f lte l ah ° nS iS reached By WW**™. ^ algorithm maps the 
patient MR images with a generic atlas without the necessity for subject-specific 
segmentation hi this way, they first have a strong confidence hT their model andmen 
progressive y take mto account the additional information coming from the data itself 
If the local refinements to be performed are small, the algorithm behaves well* 
However, if the patent's brain has been deformed too much, for example, duelo Wgh 
variability or disease the model does not work well. The algorithm works best 2 
images of normal subjects. 

t?J™fK d ™ W ™- me ^ s - de /projects/volumetry/volumetry.html, Center for 
Meorcal Diagnostic Systems and Visualisation, University of Bremen, Hahn et ai 
developed a semiautomatic extracting and volumetric analysis algorithm of the 

Z^I J™? ? fr ° m ^ image * A f6W marks were man ^ dffrneTas initial 
information for four cerebral ventricles and then the ventricular borders in 3D space 
were traced automatically by an algorithm that combined some concepts of claS 
tiansformation with a region fusion. The processing time was less than 1 second for a 
typical region of interest and the complete volumetric procedure can be performed in 
less that t minutes. Due to the complexity of the ventricle system and partial volume 
effect of MR images, the algorithm required a fine spatial resolution for the m£ 

5212?* ^ ?T 05 * addakn to the human interaction for initial 
information, the algorithm is not able to maintain the connectivity automatically. 

Lundervold A Storvik G, "Segmentation of brain parenchyma and cerebrospinal fluid 
m muluspectral magnetic resonance images," IEEE Transactions on Medical Imaging 

ITtJl t ' "t 339 " 349 deSCribeS a model - base <i segmentation method to 
extract, automatically, brain parenchyma and CSF in axial muluspectral MR images 
The algorithm incorporates information about anatomical boundaries and tissue 
signature using prior knowledge. One limitation of the algorithm is that it is restricted 
to shce images where the brain parenchyma and CSF spaces form connected regions 

DeCarli C Horwitz B, "Method for quantification of brain volume from magnetic 
ZZSl™"** U ^ Pat6nt US0052 62945A, 1993 describes a sJaZt 
T r?! 0 "^ C6ntral md ^arachnoid CSF volumes from 
MR images. The method is based on mathematical modelling of MR pixel intensity 
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histograms. The histogram is modelled as a Gaussian allowing the anoKcatiot! of 
standard statistical moments to pixel distribution. application of 

All pixels for T2-weighted MR images are considered to be CSF when the intensity of 

to for an di^ an t ValUe " ***** deViati ° n > < Wh - — 

to be set differently for different images as well as different brain structures) This is 

misleading because, m reality, one value of V or an accurate threshoTd Sep^W 

the CSF from die brain matter signals is not enough owing to the partial vo CeTffeS 

and mtensity mhomogeneity. In addition, it is difficult To locate,^ o raSuy fee 

threshold pixel mtensity from the intersection of Gaussian fitting curve? Wuse toe 

number of modals of the histogram of a ROI is difficult to determine 

wS&X ^02^f °i "* /! *** V ° lume USA patent 

US>UUo366797Bl 2002 describes a method for brain volume analysis from MRI 

images. Firstly, the brain is separated from other connected stmctoTSidZ b ^ 

surface contours are identified automatically. Then, the total SS^S me brZ 

surface contours is calculated. Fluid filled regions are thereby excluded ^d S 

accounts for the partial volume effect. The brain volume is then normafizeTby toe 

total contour volume, and toe brain parenchymal fraction is gen^ The brat 

parenchymal fraction serves as a reliable measurement of brain ftro^yld assiS to 

dete^g the severity and progression of multiple sclerosis or^conSnT 

This quantification work develops a reliable method for automatic generation oTa 

starting point for segmentation of cerebral structures, such as the laterfl ventfcL bv 

E"Auto2^ faCeS ° r d6 "5 ™ deIs > - described in Sturm B, MeieTS Hsh« 
E, Automated approximation of lateral ventricular shape in magnetic reson^ce 
images of multiple sclerosis patients," MICCAI 2002, pp P 

method it is too crude for ventricular volume analysis. however, this 

Brandt ME, "Method and apparatus for estimating tissue volumes in maenetie 
resonance ^images", USA patent US005425368A, 199? describeTa^zy app3 to 
distinguish CSF, grey matter (GM), and white matter (WM) pixeTs in m£ ont£ 
images. An unsupervised fuzzy clustering procedure based on a vLiation oftoe tozzv 
c-means algontom computes automatically toe percentage area of each of m^e ttoee 
compartments in each image. Although toe approach does not assume any a pta 

^f differed ^S^^ ^ *" ^*° f ^ 

amo^ oToverZT ^ ^ and a P^^r which determines toe 
amount of overlap of compartment boundaries. Additionally, two or more different 
spectral channels of toe same MR images are needed so £ to get more accuse 
toscrmmiation of tissue types. The approach requires about two mtouteT to analy^ a 
single brain MR unage to yield a decision on toe percent of GM, WM, a^d CSF 

Gosche KM, "Method and system for knowledge guided hyperintensitv detection m(i 
volumetric meas^ment," USA patent X^SaO^ST^S^^ 

dataset of toe brain. The system applies a vahdity-guided clustering segmentation 
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technique suitable for the discovery of small clusters to classify the varying pixel 
intensities into separate groupings, which potentially correspond to different tissue 
types. Then, the system refines the initial segmentation results into the separate 
groupings using one or more knowledge rules that combine pixel intensities with the 
spatial relationship of anatomical structures to locate one or more anatomical regions 
of the brain. Although the system uses one or more knowledge rules that combine 
pixel intensities with spatial relationship of anatomical structures to locate one or 
more anatomical regions of the brain, the system is suitable to detect brain lesion but 
not the ventricular system using multi-spectral MRI dataset. 

Despite the numerous approaches proposed to solve the problem of extracting the 
human cerebral ventricular system, the above-described methods suffer from 
disadvantages and, as such, are not suitable for clinical use. This is due to the 
necessity for human intervention, the inability of these methods to handle pathological 
and abnormal cases, and/or the inability to extract the complete ventricular system 
Moreover, the existing methods are too slow to be accepted clinically Due to the 
anatomical complexity of the cerebral ventricular system and the lack of fast and 
reliable segmentation procedures, a fast, automatic, accurate, and robust method is 
desirable to extract the complete ventricular system. 

The present invention is directed to ameliorating or overcoming the above problems 
of pnor art methods. " 

Summary of the invention 

According to a first aspect of the invention there is provided a method for extracting 
cerebral ventricular system information from images of one or more cerebral 
ventricular regions, the cerebral ventricular system comprising a third ventricle (V3) 
a fourth ventricle (V4), a left lateral ventricle (VLL) having a body (VLL-B) and an 
inferior (temporal) horn (VLL-I), a right lateral ventricle (VLR) having a body (VLR- 
B) and an inferior (temporal) horn (VLR-I), an anterior commissure (VAC) a 
posterior commissure (PC), and a midsagittal plane (MSP), the method comprising 
the steps of: r 6 

1) defining multiple regions of interest (ROI) in the images; 

2) defining seed points within each ROI; 

3) growing images of ventricular regions while correcting for leakages into 
extraventricular space; and 

4) connecting the ventricular regions grown. 

According to a second aspect of the invention there is provided an apparatus arranged 
to perform the above method. In a preferred embodiment, the apparatus is a 
computing apparatus. 

According to a third aspect of the invention there is provided a computer program 
product comprising computer program instructions readable by a computer apparatus 
to cause the computer apparatus to perform the above-defined method 
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According to a fourth aspect of the invention there is provided a method of 
quantifying the ventricular system by counting the number of voxels of the ventricular 
system extracted by a method as defined above and multiplying this count by voxel 
volume. 

According to a preferred embodiment, the ventricular system is divided into six, 
three-dimensional (3D) regions: body of the left lateral ventricle, body of the right 
lateral ventricle, inferior horn of the left lateral ventricle, inferior horn of the right 
lateral ventricle, third ventricle, and fourth ventricle. Within each ventricular region, a 
suitable plane is selected and a region of interest (ROI) is defined on it based on 
anatomy and literature studies capturing variability of the ventricular system. Within 
this ROI, its statistics are calculated which characterize the distribution of radiological 
properties (intensities) of the image data. In this invention, numerous ROIs are placed 
in some well-defined locations to study locally the distribution of cerebrospinal fluid 
(CSF), grey matter (GM), and white matter (WM). This makes the algorithm more 
robust and flexible. If the required distribution in a considered ROI is not achieved, 
this ROI is adaptively reduced and the distribution is recalculated. As the algorithm 
handles explicitly the distribution of image data, this makes it applicable to extract the 
ventricles from any type of image data. This distribution is used subsequently for 
growing each ventricular region. When a leakage into the extraventricular space 
occurs during region growing, the growing conditions are tightened in terms of 
distribution, direction of growing, and space, and growing is repeated. Finally the 
ventricular regions grown are connected by relaxing, if necessary, the growing 
conditions and applying anatomical knowledge. 

Further preferred features are defined in the dependent claims. 

Brief Description of the Drawings 

Preferred embodiments of the present invention will now be described by way of 
example and with reference to the accompanying drawings in which: 

FIG. 1 is a flowchart showing the steps in an embodiment of a method according to 
the invention; 

FIG. 2 is a radiological image showing a region of interest in a cerebral ventricle for 
calculating the radiological property distribution and seed point determination for the 
left lateral ventricle; 

FIG. 3 is a radiological image showing a region of interest in a cerebral ventricle for 
seed point determination for the inferior (temporal) horn of the right lateral ventricle; 
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FIG. 4 is a radiological image showing a region of interest in a cerebral ventricle for 
calculating the radiological property distribution and seed point determination for the 
third ventricle; 

FIG. 5 is a radiological image showing a region of interest in a cerebral ventricle for 
calculating the radiological property distribution and seed point determination for the 
fourth ventricle; 

FIG. 6a is a histogram illustrating the ideal intensity distribution of cerebrospinal fluid 
(CSF), grey matter (GM), and white matter (WM) within a region of interest; 

FIG. 6b is a part of the histogram of FIG. 6a showing the niinimum, mean, and 
maximum values of intensity for a considered component; 

FIG. 7 shows original and smoothed histograms of the intensity of matter in a region 
of interest; 

FIG. 8 is a histogram of a multiple polyline fitting of the intensity histogram of FIG. 7 
and an initial approximation of thresholds; 

FIG. 9a is a radiological image showing the division of a third ventricle (V3) into 
subregions; 



FIG 9b is a radiological image showing the division of a fourth ventricle (V4) into 
subregions; 

FIG. 10a is a schematic illustrating region connection in which first and second 
regions A and B are not connected through a third region C; 

FIG. 10b is a schematic illustrating region connection of the three regions of FIG. 10a 
showing region growing with a relaxed range; 

FIG. 10c is a schematic illustrating region connection between the regions A and B of 
FIG. 10b after they are grown towards each other following performance of the 
operations of FIG. 10b; and 

FIG. 11 is a radiological image showing how the connection between the body of a 
lateral right ventricle and inferior horn of the lateral right ventricle is achieved. 

Detailed D escription Of Preferred Embodiments 

1. Abbreviations used in the disclosure 

The following abbreviations are used in the present specification: 

- AC - anterior commissure 

- PC - posterior commissure 
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- MSP - midsagittal plane 

- VL - lateral ventricle 

- V3 - third ventricle 
V4 - fourth ventricle 

- VLL — left lateral ventricle 

- VLR - right lateral ventricle 

- VLL-B - body of the left lateral ventricle (including the anterior and posterior 
horns) 

- VLR-B - body of the right lateral ventricle (including the anterior and posterior 
horns) 

- VLL-I - inferior (temporal) horn of the left lateral ventricle 

- VLR-I - inferior (temporal) horn of the right lateral ventricle 

- AC-PC - intercommissural (axial) plane passing through the AC and PC 

- VAC - coronal plane passing through the AC 

- VPC - coronal plane passing through the PC 
CSF - cerebrospinal fluid 

- GM - grey matter 

- WM - white matter 

- ROI - two-dimensional region of interest 

- 2D - two-dimensional (or two dimensions) 

- 3D -three-dimensional (or three dimensions). 

An embodiment of the present invention provides a fast, automatic, accurate, and 
robust method to extract the complete cerebral ventricular system of, for example, a 
human being, from radiological images. The steps constituting the method are shown 
in the flow diagram of FIG, 1. The processing of the method includes at least four 
steps. The method of the present invention will be discussed in more detail after a 
brief discussion of these steps. 

The processing steps, illustrated in the flow diagram of FIG. 1, are as follows: 

1 . The first step is to load and reformat the image data. 

2. The second step is to define six 2D regions of interest (ROIs), each corresponding 
to one of the ventricular regions, VLL-B, VLR-B, VLL-I, VLR-I, V3, and V4. 

3. The third step is, for each ROI, to: 

3.1 Calculate statistics, such as the distribution of image intensities (Step 3.1). 
If for a considered ROI the required distribution (for example, two or three 
identifiable peaks in the derived histogram) is not achieved, this ROI is 
reduced adaptively and the statistics recalculated (Step 3.1a). This process is 
repeated until the required distribution is achieved. 

3.2 Define a seed point within the ROI (Step 3.2). 

3.3 Grow a 3D ventricular region from this seed point (Step 3.3). 

3.4 Check the occurrence of leakage for the ventricular region grown (Step 3.4). 
If leakage into the extraventricular space occurs, tighten the growing range in 
terms of spatial and radiological characteristics (Step 3.4a) and repeat 
growing this ventricular region. 
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4. The fourth step is to connect the regions grown: VLL-B with VLL-I, VLR-B with 
VLR-I, and V3 with V4 (Step 4). Should the corresponding regions not be 
connected, the growing range should be relaxed in terms of spatial and 
radiological characteristics and growing continued until connections are achieved 
(Step 4a). 

FIGs. 2 to 5 show radiological images of a human cerebral ventricular system. The 
ventricular system is divided into six three-dimensional regions: the body of the left 
lateral ventricle 10, as shown in the radiological image of FIGs. 2 and 3, the body of 
the right lateral ventricle 12, as shown in the radiological image of FIGs. 2 and 3, the 
inferior horn of the left lateral ventricle (not shown), the inferior horn of the right 
lateral ventricle 14, as shown in the radiological image of FIG.3, the third ventricle 
16, as shown in the radiological image of FIG. 4, and the fourth ventricle 18 as 
shown in the radiological image of FIG. 5 . 

Within each ventricular region, a plane (either axial, sagittal, or coronal) is selected 
and a region of interest (ROI) is defined on the plane based on anatomy and literature 
studies regarding variability of the ventricular system. FIGs. 2 to 5 show ROIs for the 
different ventricles illustrated, as denoted by the reference numerals 20, 22, 24 and 26 
respectively. Numerous ROIs are placed in well-defined locations to study locally the 
distribution of cerebrospinal fluid (CSF), grey matter (GM), and white matter (WM). 

According to a preferred embodiment of the invention, within an ROI, its statistics are 
calculated which characterize the distribution of radiological properties (intensities) of 
the image data. FIGs. 6a, 6b, 7 and 8 illustrate various intensity distributions of CSF 
GM and WM within an ROI. These are discussed in more detail below. If the 
required distribution (two or three identifiable peaks in the intensity histogram) in a 
considered ROI is not achieved, this ROI is adaptively reduced and the distribution is 
recalculated. This distribution is used subsequently for growing each ventricular 
region. When a leakage into the extraventricular space occurs during region growing, 
the growing conditions are tightened in terms of distribution, direction of growing' 
and space, and growing is repeated. Finally the ventricular regions grown are 
connected by relaxing, if necessary, the growing conditions and applying anatomical 
knowledge. 

To facilitate region growing, the ventricular regions may be further subdivided into 
smaller subregions, as shown in FIGs. 9a and 9b and denoted by reference numerals 
1,2, 3 and 4. FIG. 9a shows the division of a third ventricle (V3) into subregions. V3 
is divided into four subregions by the planes passing through AC and PC, namely 
VAC, VPC, AC-PC. The subregion 1 contains the most antero-ventral part of V3 
separated from the rest by VAC and AC-PC planes. The subregion 2 contains the 
most dorsal part of V3, separated from the rest by AC-PC and VPC planes. The 
subregion 3 contains the most posterior part of V3, separated from the rest by AC-PC 
and VPC planes. The subregion 4 contains the most ventral part of V3, separated from 
the rest by AC-PC and VAC planes. 
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FIG 9b shows the division of a fourth ventricle (V4) into subregions. V4 is divided 
into two subregions by an axial plane 30 passing through the point from which the 
growing takes place, namely the seed point. The suhregion 1 includes the superior 
part of V4, and the suhregion 2 includes the inferior part of V4. 

FIG. 10a shows region connection in which first and second regions A and B are not 
connected through a third region C. FIG 10b illustrates region connection of the 
three regions of FIG. 10a. Local statistics are calculated in rectangular areas 32 at 
both ends of regions A and B, and region growing with a relaxed range is performed. 
FIG. 10c illustrates region connection in which both regions A and B are grown 
towards each other after performing operations specified in FIG 10b, but regions A 
and B are still not connected. New rectangular areas 34 are set, statistics calculated 
within them, and the region growing is performed again. 

FIG. 11 illustrates the connection between the body of a lateral right ventricle and 
inferior horn of the lateral right ventricle which is discussed in detail below. 

2. Algorithm 

2 -l STEP 1 :Loading and reformatting of image tlata 

The image data of the brain are loaded first. It is assumed that the location of the 
MSP and the positions of the AC and PC are known. The location of the MSP may be 
calculated, for example, using the method described in International Patent 
Application No: PCT/SG02/00006 filed Jan. 2002 in the names of Hu Q, and 
Nowinski WL, and entitled "Method and apparatus for determining symmetry in 2D 
and 3D images". The positions of the AC and PC may be calculated, for example, 
using the method described in International Patent Application No: PCT/SGOO/00185 
filed Nov. 2000 in the names of Nowinski WL, and Thirunaruukarasuu A, and entitled 
"Methods and apparatus for processing medical images". 

As the MSP, the AC and the PC are known, the image data of the brain are 
reformatted along the AC-PC plane, the MSP or the VAC such that axial, coronal, and 
sagittal orientations are available (Step 1 shown in FIG. 1). This permits the use of 
any orientation for setting an ROI and performing oriented region growing. 

If the image data are acquired with low inter-slice resolution, for example, 5 mm or 
more, reformatting does little to assist in achieving a better accuracy of segmentation 
and quantification. In such an instance, reformatting should be carried out only to 
determine ROIs and several of the refinements and tests proposed below may not be 
applicable. However, region growing for all regions should be attempted in the 
acquisition orientation. 

The low resolution case is an unusual situation and the accuracy of obtaining 
quantification of the ventricular system from low resolution data may be unreliable. 



WO 2004/077359 



PCT/SG2003/000043 



13 



2.2 STEP 2 : Defin e multiple regions of interests (ROIs) 

mTTmo?^ C* 018 ) 316 de ^d corresponding to the ventricular regions: 

3^3 . V ^" B ' VLL_I ' VLR - 15 V3 ' ^ V4 > < Ste P 2 > FIG - D- ^ initial ROfa are 
denned taking mto account the worst case assumption, that is, the biggest ROI 
necessary, based on literature studies. " 

2.2.1 ROIs for VLL-B and VT/R-R 

The ROIs are set on the coronal orientation as follows 

1. Select the VAC. 

2. Set the initial rectangular ROI on the VAC: 

f) l^™ 11 / between [MSP, MSP + 25 mm] for VLL-B and [MSP-25mm, 
MSP] for VLR-B, and 

ii) dorsally between [AC, AC + 35 mm] (the sign "+" means in upper 
(dorsal) direction from the AC). 

3. Refer to 2.3.2 below for the adaptive determination of ROI for subsequent 
processing. 

Both lateral ventricles (10 and 12 as shown in FIGs. 2 and 3) are very variable in 
shape and volume. The average volume is 7 ml but the volume may vary from around 
2 to 39 ml. The estimates of the size of lateral ventricles 10 and 12 given in prior art 
literature are based on various techniques such as ventricular casts and 
pneumoencephalograms which may not be accurate. The left lateral ventricle 10 is 
slightly larger than the right ventricle 12 in most cases. The distance between the 
lateral margins of both ventricles at its minimum spread is typically 30 to 40 mm (see 
Radiology of the skull and brain, ventricles and cisterns. Eds: Newton TH and Potts 
DG. MediBooks Great Neck, NY, pp. 3494 - 3537). Thus, the rxiinimal distance 
between the MSP and the lateral edge of either of the lateral ventricles is around 15 to 
20 mm, so the lateral ROI range of MSP to MSP + 25 mentioned above would appear 
to be justified as it encloses the lateral boundary of the lateral ventricles. 

Although the height of the ROI is not generally available in the prior art literature it 
can be calculated from the method of Jirout [see Robertson EG 
Pneumoencephalography, ed. 2, Springfield H. 1967, Thomas CC, Publisher] The 
height of the lateral ventricles can be determined and is typically around 32 to 34 mm 
(if angle (the angle between the tangents drawn along the superomedial and 
mferolateral walls of lateral ventricles in coronal plane) is in the range 20° to 25° ) 
Thus, a height of AC+35 would also appear to be justified. 

2.2.2 ROIs for VLL-I and VLR-I 

The temporal horns are very variable. The maximum spread of the terminal (most 
ventral) points is around 57 mm and the width of the temporal horns is around 15 to 
20 mm They curve medially as they progress posteriorly and the maximum spread is 
around 73 to 90 mm at the maximum convexity [see Radiology of the skull and brain 
ventacles and cisterns. Eds: Newton TH and Potts DG. MediBooks, Great Neck NY 
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^foUows: 20 ^ 22> 38 mUStKlted * FIGS - 2 and 3 « Set ° n 1116 coronal orientation 
1 . Select all the coronal slices between the VAC and the VPC 

2 * ]?SiSL? ?P 7£ ° f the 4 anterior P 3 * 8 of *• tem P°ral horns to that of the widlh 
££^£ • 5 l aCC ^ dm 8 to Radiology of the skull and brain, ventricles and 
cisterns. Eds: Newton TH and Potts DG. MediBooks, Great Neck, NY, pp. 3494 

m 11' * ? W 6 1 h 1 onzontal P rofile at the level of the inferior point of V3. 
Identify the temporal bone signal in the profile as the lateral-most high signal due 
to the bone marrow fat and measure the distance from this to the MSP. Divide 

25J£ TTZ m 5^ ° 0nSider MS 38 P ° mt " a " ™ s P oint should ^ in the 
vicinity of the ventral part of the temporal horn. Draw the ROI 20 mm x 20 mm 

in reference to this point so that the coordinates are (a, a-10) fa, a+10^ laterallv 
and (a, a-20) ventrally (inferiorly), [see FIG. 3] ^ * V 

3. Keep unchanged the ROI for all subsequent processed slices. 

2.2.3 ROI for V3 

Stow? 1 ^ V3 ' ™ illUStrated fa FIG - 4 > is set on orientation as 

1. Select the MSP. 

2. Set the initial rectangular ROI: antero-posteriorly between [AC, PC] and dorsally 
between [AC -10 mm, AC + 35 mm], [see FIG. 4]. X 

3 " proSss^g 2 ' 3 ' 2 b6lOW ^ ada P 6vc determination of ROI for subsequent 

2.2.4 ROI for V4 

Initially a large ROI 26 is used for V4, as illustrated in FIG. 5. The length of 
aqueduct is token to be around 12 to 18 mm, and the length of floor of V4 as around 
29 mm and its height as around 20 mm maximum according to Radiology of the skull 

SSFSy V T£Z "J^fTi EdS: TH md P0ttS DG - Medfflooks, Great 

W ^' 0 PP ' ? " 1 3537, • V 18 SUmlar ^ sha P e to a rid S ed te *> that is, it has a 
base having a rectangular penpheral outline and sloping sides extending therefrom to 
an apex pointing at the fastigium posterosuperiorly. This is shown in FIG. 5. This 
ROI may include , the clivus (bone), the ambient cistern, quadrigeminal cistern, and 
cistema magna. These undesirable structures can be excluded if the ROI is made 
triangular and an apex pointing at the fastigium posterosuperiorly. 

The ROI is set on the sagittal orientation as follows- 

1. Select the MSP. 

2. Set the rectangular region: 

i) posteriorly between [PC, PC + 50 mm] (the sign '+' means in posterior 
direction with respect to the PC), and 

ii) dorsally between [PC, PC - 45 mm] (the sign means in ventral 
direction with respect to the PC). 
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3. Limit the superoposterior part of the ROI by drawing a line starting at 5 mm 
posterior to the PC to the posteroinferior corner of the ROI to bisect the ROI 
obliquely. 

4. Limit the inferoposterior part of the ROI by drawing the line from the 
inferoantenor comer of the ROI to the midpoint of the posterior limit of the ROI. 
This will change the original square ROI into a pentagon with a base facing 
anteriorly and its apex roughly corresponding to the apex of V4 [see FIG. 5]. 

5. Refer to 2.3.2 below for the adaptive determination of ROI for subsequent 
processing. 

2.3 STEP 3 - Calculate st atistics for each ROT 

The calculated statistics characterize the distribution of intensities in the studied 
region of interest. The initial ROIs 20, 22, 24 and 26 are chosen such that each 
contains three components: CSF, GM, and WM, (Step 3.1, Figure 1). For a considered 
ROI, its histogram is calculated and the gray level ranges of CSF, GM, and WM are 
determined, as illustrated in FIG. 6a. For each component, three values are calculated- 
mean, minimum and maximum, as shown in FIG. 6b. The minimum and maximum 
values correspond to valleys while the mean values correspond to the peaks in the 
histogram. Histograms of real data have substantial local variations, which hampers 
determination of global peaks and valleys. Therefore, some initial smoothing can be 
apphed and intensity distribution approximated, for instance, by employing 1) 
triangular, 2) polyline or 3) Gaussian approximations. 

2.3.1 Histogram processing 

Each ROI 20, 22, 24 and 26 is defined to facilitate analysis that will produce 
statistically significant findings necessary for ventricular system extraction. FIG. 7 
shows the histogram obtained from an ROI containing the lateral ventricles and 
adjacent structures. This is designed to produce a multi-modal histogram where the 
intensity peaks can be chosen, corresponding to WM, GM, CSF (and other classes if 
necessary). The ventricles are defined by those spatially connected pixels with 'an 
intensity between two thresholds, so called intensity thresholds, as described below. 

There are five steps to determine the intensity thresholds from histogram peaks and 
valleys: 

1) calculating a specific histogram of ROI; 

2) smoothing the histogram; 

3) identifying the peaks in the histogram; 

4) determining the peaks corresponding to intensity values of CSF, GM, and WM 
using neuroanatomical knowledge; and 

5) determining the thresholds. 

FIG. 7 shows the original and smoothed histograms. The smoothed histogram is 
obtained by applying Fourier Transformation to the original histogram. The 
following constraints are imposed to exclude undesirable peaks: 
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1. A peak is ignored if its area is small compared to the cumulative area of the 
histogram (i.e., if its area is less than, for example, 1730 th of the total area) 

2. A peak is ignored if it is between the intensities at, for example, 5% and 95% 
of the cumulative area of the histogram. 

After finding all of the peaks in the histogram, neuroanatomical and radiological 
knowledge is incorporated to locate the peaks corresponding to CSF, GM and WM 
(for example, it is assumed that the highest intensity peak corresponds to WM on Tl- 
weighted MRI while the lowest intensity peak corresponds to CSF). 

Most existing methods calculate the CSF-GM threshold by taking the average of the 
two peaks corresponding to CSF and GM, and calculate GM-WM thresholds by 
taking the average intensity of the peaks corresponding to GM and WM. In this 
invention, the initial thresholds are approximated using polyline or Gaussian fitting 
The initial thresholds are approximated by polyline fitting of the histogram illustrated 
as follows: 

• Smooth the histogram of ROI, as shown by line 1 in FIG. 8. 

® In the smoothed histogram, locate the positions of peaks and valleys, and 
connect the peaks and valleys sequentially, as shown by line 2 in FIG 8 
(first polyline fitting). 

• For each line segment of the first polyline fitting, locate the furthest point 
from the smoothed histogram to this line segment. 

• Connect all the peaks, valleys, and the furthest points of the smoothed 
histogram, as shown by line 3 in FIG. 8 (second polyline fitting). 

© Locate the peaks or intensity values corresponding to CSF and WM guided 
by neuroanatomical and radiological knowledge (i.e., it is assumed that the 
highest intensity peak is of WM on Tl-weighted MRI. The lowest 
intensity peak is of CSF). 

• Locate the previous furthest point and the next furthest point for the peaks 
of CSF and WM on the second polyline fitting. These points are the initial 
thresholds for CSF-GM, GM-WM, and so on, and correspond to the 
vertical lines shown in FIG. 8. 

These initial thresholds will be changed adaptively during the growing of the 
ventricular system to cope with the partial volume effect and intensity inhomogeneity. 

2.3.2. Reduce ROI adaptively 

As the initial ROIs 20, 22, 24 and 26 are defined based on the worst case prior art 
literature variability studies, the distribution of CSF, GM, and WM may be 
unbalanced, in some cases, (meaning that some peaks corresponding to either WM 
GM or CSF are so small that they are not easily identifiable in the histogram) This 
may hinder an accurate processing of the histogram. To improve the density 
distribution balance, the ROI is reduced adaptively. When the ratio of the WM/GM to 
CSF peak is greater than a given number (4, for example), the ROI is reduced by a 
given percentage (20%, for example), and the histogram is recalculated. The adaptive 
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reduction of ROI is performed iteratively until the required distribution is achieved 
(Step 3.1a, Figure 1). 

For the ROIs corresponding to VLL-B and VLR-B, the aspect ratio is kept constant, 
and the ROIs are reduced from the lateral and dorsal sides. The ROI 24 corresponding 
to V3 is reduced from the superior side. The ROI 26 corresponding to V4 is reduced 
from the posterior side only. The ROIs corresponding to VLL-I and VLR-I are not 
used for statistics calculations and consequently remain unchanged (see 2.4.2 below). 

2,4 Define seed point within ROI 

A single seed point is defined within each ROI (Step 3.2, Figure 1) and the region 
growing starts from this seed point. 

2.4.1 Seed points for VLL-B and VLR-B 

Each of the seed points for VLL-B and VLR-B is defined as follows. 

1 . Sample the intensity of the ROI horizontally, for example every 5mm starting 
fromAC+lOmm. 

2. Calculate a profile along each sample line segment. To increase the robustness, 
the averaged profile with, for example, two additional lines passing one pixel 
below and one pixel above is used. 

3. Find the longest CSF segment from the averaged profile and place the seed point, 
for example, in the middle of this segment 

2.4.2 Seed points for VLL-I and VLR-I 

The exact position of the right coronal plane containing the ROI with respect to the 
reference planes (VAC, VPC, and AC-PC plane) is unknown. Therefore, the ROI 
defined in 2.2.2 above is movable and it is set on multiple coronal planes between the 
VAC and the VPC. Planes are selected and processed starting from the VAC 
posteriorly. The CSF region on the coronal slice between the VAC and the VPC with 
the largest CSF region is located. As the statistics are not known in this ROI, the 
VLL-B statistics with the CSF range are used for VLL-I and that of VLR-B for VLR- 
L The seed point is placed at the geometric center of this largest CSF region. 

2.4.3 Seed point for V3 

The seed point for V3 is defined as follows: 

1) Calculate a profile along the AC-PC line. 

2) Determine a pixel on the AC-PC line segment, whose gray level is closest to the 
CSF mean value, and take it as the seed point. 

2.4.4 Seed point for V4 

The seed point for V4 is defined as follows: 

1 . Sample the intensity of the ROI horizontally, every 2 mm, as described in 2.2.4 
above, starting for example, from PC- 1 0 mm. 
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2. Calculate a profile along each sample line segment (to increase the robustness 
several lines, 3 for example, can be used to obtain the averaged profile as in 2 4 1 
above). 

3. Calculate the length of the CSF in each averaged profile and compare the length to 
the previous one. When this length starts decreasing for at least two subsequent 
fine segments, take, for example, the middle of the longest CSF segment as the 
seedpomt. 

2.5 Grow each ventricular region 

The ventricular regions are grown in 3D independently starting from the defined seed 
points (Step 3.3, Figure 1). Region growing is directional which allows for better 
control of growing in 3D space. 

Let m be the minimum, M the maximum and the mean values of the CSF range 
calculated m Step 3.1. By using the complete range of intensities [m, M], the region 
grown may be overestimated because of the partial volume effect. Let s be a scaling 
factor between 0 and 1. Region growing can then be better controlled by using the 
foUowmg growing range [ -s*( - m ), +s*(M- )] with a variable value of s. For 
J," 0 ! onl y mean val ue of CSF is used for growing. For s = 1, the full range of 
CSF is utilized. For s = 0, the region grown may be underestimated while for s=l it 
may be overestimated. The value of s has to be selected based on quantitative 
assessment. 

To facilitate region growing, the ventricular regions are further subdivided into 
smaller subregions, as illustrated in FIGs. 9a and 9b. This approach has several 
advantages, namely: 

• Region growing is simplified as complex shapes are replaced by simpler ones. 

• Easier control regarding growing and connecting. 

• Better leakage control, as it is easier to incorporate specific domain knowledge 
m each subregion. 

• Processing is more efficient as only a subregion needs to be regrown in case of 
leakage. 

• Facilitated reduction of the partial volume effect, as it is easier to incorporate 
specific domain knowledge in each subregion. 

• Easier to adjust the initial thresholds tailored to the local anatomy. 

2.5.1 Growing of VLL-B and VLR-B 

Each of the VLL-B and the VLR-B regions is grown in 3D space on coronal slices 
shce by slice. Growing is initiated anteriorly from the seed point located on the VAC 
When this growing is completed, it is continued posteriorly on all subsequent coronal 
slices. Eventually, it is continued anteriorly when attempting to extract the posterior 
part of the inferior horn. 
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2.5.2 Growing of VLL-I and VLR-I 

Each of the VLL-I and the VLR-I regions is grown in 3D space on coronal slices, 
slice by slice. Growing is initiated anteriority from the seed point. When this is 
completed, growing is continued posteriorly on subsequent coronal slices. 

2.5.3 Growing ofV3 

V3 is divided into four subregions by the planes passing through the AC and the PC, 
namely the VAC, the VPC, and the AC-PC, as shown in FIG.9a. The subregion 1 
contains the most antero-ventral part of V3, separated from the rest by the VAC and 
the AC-PC planes. The subregion 2 contains the most dorsal part of V3, separated 
from the rest by the AC-PC and the VPC planes. The subregion 3 contains the most 
posterior part of V3, separated from the rest by the AC-PC and the VPC planes. The 
subregion 4 contains the most ventral part of V3, separated from the rest by the AC- 
PC and the VAC planes. 

V3 is grown in 3D as follows. Subregion 2 is grown dorsally and subregion 4 is 
grown ventrally on axial slices starting from the seed point. Subregion 1 is grown 
anteriorly on coronal slices starting from any pixel common with subregion 4. 
Subregion 3 is grown dorsally on axial slices starting from any pixel common with 
subregion 2. 

2.5.4 Growing V4 

V4 is subdivided into two subregions by the axial plane passing through the seed 
point, as shown in FIG. 9b. Subregion 1 includes the superior part of V4 and 
subregion 2 includes the inferior part of V4. 

V4 is grown on axial slices, dorsally in subregion 1 and ventrally in subregion 2 
starting from the axial slice containing the seed point. 

2.6 Leakage conditions for each region grown 

The foreground region is formed by all pixels belonging to the extracted ventricular 
system. The remaining pixels, i.e., those not belonging to the ventricular system, form 
the background region. 

First the leakage has to be detected. Leakage in various areas behaves in different 
ways so a single algorithm does not work for all the regions. Anatomically there is no 
sudden change in the size of the ventricular system except in specified places. 
Therefore, when a region is grown on a given slice it is compared with the previous 
slice and any increase in size more than, for example 50%, can be considered a 
leakage. If the width of V3 is, for example, 4 mm in one slice and if it exceeds, for 
example, 6 mm (50%) in the subsequent slice then leakage is considered to have 
occurred. At other locations, a grown region can be considered leakage if it is grown 
to an area that is not anatomically correct. 
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2*6.1 Anti-leakage conditions 

A region grown may "leak" beyond the ventricular space into the extraventricular 
space. Two groups of conditions are formulated to deal with this problem: 1) leakage 
detection conditions and 2) leakage control conditions. These conditions are based on 
domain knowledge (anatomy and variability studies), radiological properties 
(intensity ranges), geometrical constraints, direction of growing, and size and location 
of the grown region. 

2.6.2 VLL-B and VLR-B leakages 

Leakage may occur in a posterior ventral part of the body of lateral ventricles. 
Moreover, in some cases, a large cystic elongated space between the lateral ventricles 
may be present representing the cavum septum pellucidum and cavum vergae. A 
potential leakage to this space has to be checked as well. This can be checked on 
radiological images as the cavum is relatively perpendicular and predictable, and the 
septal wall can be easily detected (it has higher signal intensity as compared to CSF 
on Tl-weighted images). For the leakage into the quadrigeminal cistern, any growth 
of a region below the PC, in the area behind the VPC up to the seed point of V4, can 
be considered as leakage. 

2.6.3 VLL-I and VLR-I leakages 

The following leakages may occur in these regions: 1) inferomedially through the 
hippocampus to the interpeduncular cistern, 2) medially through white matter (optic 
radiation) and cortex (choroidal fissure) to the subarachnoid space. Suitable spatial 
constraints are imposed for leakage prevention. This is described in more detail in 
2.6.6.2 below. 

2.6.4 V3 leakages 

The following leakages may occur in V3: 1) anteriorly through the lamina terminalis 
to the chiasmatic cistern, 2) ventrally through the mesencephalic tegmentum to the 
interpeduncular cistern, 3) posteriorly through the posterior commissure (stalk of the 
pineal body) to the cisterna ambiens. A leakage is considered to have occurred when 
the number of extracted pixels in the current slice is, for example, 50% more than that 
in the previous slice. 

2.6.5 V4 leakages 

The following leakages may occur in V4: 1) dorsoposteriorly through the superior 
medullary velum to the cisterna ambiens, 2) ventroposteriorly through the inferior 
medullary velum to the cisterna magna. A leakage is considered to have occurred 
when the number of extracted pixels in the current slice is, for example, 50% more 
than that in the previous slice. 

2.6.6 Tighten growing range 

When the leakage occurs, the region growing has to be repeated with tightened 
growing ranges. The range can be narrowed from intensity and spatial points of view. 
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By decreasing the value of the scaling factor s, the growing range of intensity [ -s* 
( -m), +s*(M- )] is narrowed. The tightening is done iteratively (Step 3.4a) by 
decreasing the value of s (for example by 5% starting from s = 0.9), followed by 
performing the region growing again and checking whether the leakage is avoided. 

When the increase in the number of extracted pixels from a slice exceeds the number 
from the previous slice by, for example, more than 50%, it can be considered that 
leakage has occurred. When this happens, growing is tightened iteratively by 
decreasing the scaling factor s. This process is repeated till the increase in the number 
of extracted pixels in a slice does not exceed that of the previous slice by, for 
example, more than 50%. Alternatively, growing may be tightened iteratively by 
decreasing the CSF-GM threshold ranges. This process is repeated till the increase in 
number of extracted pixels is less than, for example, 50% in the subsequent slice. 
Spatially, if leakage has occurred, a bounding box is decreased as follows to confine 
the growing. 

2.6.6.1 VLL-B and VLR-B 

Leakage may occur in the quadrigeminal cistern of the body of lateral ventricles 
behind the VPC up to the seed point of V4. The growing will not be allowed below 
the PC in ventral direction. 

2.6.6.2 VLL-I and VLR-I 

To prevent potential leakages, the foreground region(s) on a processed coronal slice is 
(are) not allowed to grow beyond the boundary of ROI for VLL-I or VLR-I. 

2.6.6.3V3 

As the shape of the roof of V3 is oblique, the distance between the superior and 
inferior pixels of the foreground region decreases along the anterior direction. In order 
to prevent the potential leakage ventrally through the mesencephalic tegmentum to the 
interpeduncular cistern, the width of the foreground region of subregion 1 on a 
processed coronal slice is not allowed to increase by, for example, more than 50% of 
the width of the previous slice along the anterior direction. 

To prevent leakage posteriorly through the PC (stalk of the pineal body) to the 
cisterna ambiens, two spatial constraints are imposed: 

• The maximum width of the foreground region of subregion 3 on a currently 
processed axial slice is on the PC line (the PC line is the intersection between 
the VPC and the current axial slice). 

• The distance between the gravity centre of the foreground region of 
subregion 3 and the MSP should be small, for instance less than 4 mm [see 
Radiology of the skull and brain, ventricles and cisterns. Eds: Newton TH 
and Potts DG. MediBooks, Great Neck, NY, pp. 3494 - 3537]. 
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2.6.6.4 V4 

To prevent leakage dorsoposteriorly through the superior medullary velum to the 
cisterna ambiens, two constraints are imposed: 

• The number of the foreground pixels of the current axial slice should not 
increase by, for example, more than 50% over that in the previous slice in the 
inferior direction. 

• The distance between the gravity centre of the foreground region and the MSP 
should be less than, for example, 2 mm. The gravity centre of the foreground 
region should not deviate, for example, more than 2 mm away from that of the 
previous slice. 

To prevent the leakage ventroposteriorly through the inferior medullary velum to the 
cisterna magna, the width of the foreground region at the lateral recesses should be 
less than, for example, 2 mm [see Radiology of the skull and brain, ventricles and 
cisterns. Eds: Newton TH and Potts DG. MediBooks, Great Neck, NY, pp. 3494 - 
3537J • 

2.7 STEP 4 - Conn ect the corresponding grown regions 

The ventricular regions grown, VLL-B, VLR-B, VLL-I, VLR-I, V3, and V4 may 
generally, not be connected as the CSF regions between them are quite narrow and the 
partial volume effect, noise, and low image resolution may prohibit growing through 



The algorithm first checks whether the corresponding ventricular regions are 
connected. If not, the growing conditions are relaxed, and region growing is 
attempted, see as an example, FIGs. 10a, 10b and 10c. These steps are repeated to 
achieve the desired connection. If the regions cannot be connected through relaxing 
the growing conditions (i.e., based on radiological properties), (fine) anatomical 
knowledge is incorporated to ensure connection. If the regions still cannot be 
connected, brute-force conditions (based on crude anatomical knowledge) are applied. 

2.7.1 Growing region under relaxed conditions 

The gray level ranges are relaxed by increasing the value of s (for example by 5% 
starting from s = 0.5) in the intensity range determined by the formula [ -s*( -m) 
+s*(M- )] . This concept is illustrated in FIGs. 10a, 10b and 10c. To deal with 
the partial volume effect, the value of s can be greater than 1 . 

Region growing under relaxed conditions is applied to the following situations: 

• Connect the extracted regions of the ventricular system, 

• Segment regions with intensity inhomogeneity, i.e., subregion 2 of V4, 

• Extract the narrow parts of the lateral ventricle, i.e., the posterior horn.' 
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2.7.2 Connect VLL-B with VLL-I and VLR-B with VLR-I 

The detailed steps to ensure connection between VLR-B and VLR-I are set out below. 
This connection is performed in the coronal direction. The connection between VLL- 
B and VLL-I is achieved in a similar way. 

When growing in the posterior direction, the following steps are involved: 

1. Locate the leftmost inferior pixel Pi of the VLR-B. 

2. Grow VLR-I on a current coronal slice in the posterior direction. The current 
coronal slice is denoted . Check whether the number of foreground regions 
is equal to or less than 1 . If not, go to * 3 ' below. If yes, 

• Find the minimum rectangular area containing all the foreground 
pixels of the previous coronal slice in the anterior direction. 

• Copy the rectangular area to the current coronal slice. 

• Find a seed point for the grey matter at the rectangular area and extract 
the grey matter through region growing. If the seed point cannot be 
found or if the grey matter is absent, go to step '6' below for a brute- 
force connection. 

• Trace the leftmost superior boundary of the grey matter. 

• Add the traced pixels to the grown region, and go to step '4' below. 

3. Check whether the number of foreground region(s) is more than 1 on . If 

not, go to step '4' below. If yes, connect the foreground regions by" the 
following steps: 

• Locate the P 2 and P 3 pixels on the foreground regions, where P 2 is the 
leftmost inferior pixel of the upper foreground region and P 3 is the 
leftmost superior pixel of the lower foreground region, as shown in 
FIG. 11. 

• Trace the boundary of the grey matter between P 2 and P 3 . 

• Add the traced pixels to the grown region. 

4. Check whether the current grown region is connected with the grown region 
on the previous coronal slice in anterior direction. If not, connect them. 

5. Check whether Pi is on the current coronal slice. If not, go to step '2' above. If 
yes, terminate the current connection process. 

6. The brute-force connection is enforced when either the seed point is not found 
or the grey matter is absent. Form a fan-type layer from Pi and the last grown 
region to complete the connection. 

2.7.3 Connect V3 and V4 - Extraction of Aqueduct 

In most cases, the aqueduct, that is, the connection structure between V3 and V4 
cannot be extracted directly as its diameter is quite small (about 1.2 mm). After 
finishing segmentation of V3 and subregion 1 of V4, the extraction of the aqueduct is 
done in the superior direction as follows: 

1. Locate the superior aqueduct starting pixel (Qi) at the most posteroinferior 
part of V3. Locate the most superior pixel of subregion 1 of V4 (Qj), FIG 10b. 
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The axial slice with Q 2 is called reference slice, A R . The following 
assumptions are made when locating the superior aqueduct start pixel (Qi): 

• Its location is inferior to PC. 

• Its location is at the most postero- ventral part of V3 . 

• It is located near the MSP. 

2. Check whether Qi and Q 2 are spatially connected or the same. If yes, V3 is 
connected with V4, and the current connection is terminated. 

3. Determine a rectangular area around a pixel (described below) on the axial 
slice next to 4, in superior direction. 

• The intersection of the line connecting Qi and Q 2 with the axial plane 
in the superior direction is called the intersection pixel (IP). 

• The half height of the rectangle in the anterior-posterior direction 
equals to the distance between IP and Q 2 . 

• The half width of the rectangle in the left-right direction is equal to the 
distance between Q 2 and the geometric centre of the foreground region 
on the axial slice next to A R in inferior direction. 

4. Grow region on the axial slice next to A R in superior direction. The aqueduct 
can be positioned at the lowest intensity in the search area on Tl -weighted 
MR image (or the highest on T2-weighted images). If the number of pixels 
grown is zero, go to step 7 to force the connection. 

5. Check whether the grown region connects with subregion 1 of V4. If not, 
connect the grown region with the subregion 1 of V4. 

6. Add the grown region into subregion 1 of V4. Assign the current axial slice as 
the reference slice A R . Locate the IP of grown region as Q 2 , go to step 2. 

7. For the brute-force connection, take the IP in each axial slice as aqueduct. 

Ensuring connections (if they do not exist) between the lateral ventricles with the 
third ventricle (through the foramina of Monroe) may not be necessary, as these 
foramina open temporarily under the increasing pressure of CSF and close again. 
Therefore, the algorithm does not provide any additional means to ensure these 
connections automatically. The decision on ensuring these connections can be made 
interactively by the user of the algorithm, and the procedure is conceptually similar to 
that for providing the connection between V3 and V4. 



In summary, an embodiment of the present invention is directed to a method and 
apparatus for extracting the human cerebral ventricular system using a combination of 
neuroanatomy, radiological properties, variability of the ventricular system taken 
from literature, and image processing techniques used for extraction. The 
embodiments of the invention enable any type of image data to be processed and 
provide a rapid, robust, and accurate means for the extraction of the complete 
ventricular system using a fully automatic method. 
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Various modifications to the methods described above will be apparent to those 
skilled m the art, for example, whilst the present invention has been described above 
in connection with the extraction of the cerebral ventricular system for a human being 
it will be appreciated that the methods could be applicable to animals, as animal 
brains also have a similar morphology and the ventricular system is even 
comparatively larger in lower animals. 



